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A nonperturbative electron transfer rate theory is developed based on the reduced density matrix 
dynamics, which can be evaluated readily for the Debye solvent model without further approxima- 
tion. Not only does it recover for reaction rates the celebrated Marcus' inversion and Kramers' 
turnover behaviors, the present theory also predicts for reaction thermodynamics, such as equilib- 
rium Gibbs free-energy and entropy, some interesting solvent-dependent features that are calling 
for experimental verification. Moreover, a continued fraction Green's function formalism is also 
constructed, which can be used together with Dyson equation technique, for efficient evaluation of 
nonperturbative reduced density matrix dynamics. 



I. INTRODUCTION 

Solvent environment plays a crucial role in determin- 
ing chemical kinetic properties. Its interaction induces 
fluctuation that compensates for the activation energy, 
and also results in relaxation that stabilizes the reac- 
tion product. This dual role of solvent interaction, re- 
lated via fluctuation-dissipation relation, has been recog- 
nized since Einstein in his study of Brownian motion.^ 
The effect of solvent interaction on chemical kinetics was 
first studied by Kramers in his classical Fokker- Planck- 
equation (FPE) approach to the rate theory of isomer- 
ization reactioniS This is a classical reduced equation-of- 
motion (EOM) approach, in which chemical reaction is 
treated as escape over barrier of particle moving in one- 
dimensional potential surface and subject to a Langevin 
force originated from stochastic solvent interaction. The 
resulting rate, as function of solvent viscosity, exhibits a 
turnover behavior: it increases linearly as viscosity ini- 
tially, reaches a maximum at an intermediate viscosity 
value, and falls off inversely as viscosity in the high fric- 
tion regimeiSi^ 

Electron transfer (ET) processes constitute another 
class of systems whose dependence on solvent environ- 
ment has been extensively studied since Marcus' pioneer- 
ing contributioniii^ Consider the simplest ET system in 
which the total Hamiltonian reads 

H^ = ha\a){a\ + {hb + E°)\b){b\ + V{\a){b\ + \b){a\). (1) 

Here, ha and hb are the solvent Hamiltonians for the ET 
system in the donor and acceptor states, respectively, E° 
is the reaction endothermicity, and V the transfer cou- 
pling matrix element. The system is initially in the donor 
I a) site, with the solvent (bath) equilibrium density ma- 
trix p°^ ex e^'*"/^''^-'") at the temperature T. The reaction 
coordinate is now [/ = hh — ha, which is purely of solvent 
in nature and called the solvation coordinate. Its static 
mean, A — (U) = tr(C//9®'^), denotes the solvent reorgani- 
zation energy, while its variation satisfies (C/^) — (C/)^ w 



2XkBT, the classical fiuctuation-dissipation relation in 
the high temperature limit. With a classical dielectric 
continuum model, Marcus had further related the solvent 
reorganization energy A to the static and high-frequency 
dielectric constants of the solvent 4*^ 

The standard approach to ET rates is based on correla- 
tion function formalism,4i5i6,7,8,9,10,ll,12,13,14,15,16,17,18,19 

The second-order transfer correlation function is re- 
lated to the nonadiabatic rate theory,/"-^-^ which in the 
static solvation limit assumes the celebrated Marcus' rate 
expression (^i^ 



k = 



vyh 



y/XkBT/n 



exp 



{E° + A)2 



UknT 



(2) 



Rates have also been formulated based on fourth- 
order transfer correlation functions, followed by 
certain resummation schemes that partially ac- 
count for the effect of nonperturbative transfer 
couplingi7ii?iiiii^ii?ii4ii^ii?ii7ii?ii? The resulting rates do 
recover the celebrated Kramers' turnover behavior jS*^ 
and are also able to relate the reaction adiabatic- 
ity to solvent relaxation time scale. However, these 
correlation function-based rate theories remain per- 
turbative in nature. To obtain rate expressions, one 
shall also assume the separation of time scales be- 
tween the fast ET dynamics and the slow solvation 
processesj ^i^°i^^i^^'^^i-^^i^^i-^^'-^^i^^i-^^ As results, the reac- 
tion Gibbs free-energy AG", obtained via the forward 
and backward rates ratio k/k' = exp[— AG°/(A:bT)], 
is identical to the endothermicity parameter E° that 
contains no explicit dependence on solvent environment. 
The parabolic dependence of log k on E° , as shown by 
eq|21 can be read as its dependence on AG°. This is the 
so-called Marcus' inversion behavioriii^ 

Alternative approach to ET rates is via reduced density 
matrix, ?2i^ defined formally as p{t) = trBPT(t), the trace 
of total density matrix over bath degrees of freedom. This 
is a quantum reduced EOM approach in which the trans- 



fer coupling is part of the system and can be treated ex- 
actly. However, the system-bath interaction, which in 
ET systems is neither weak nor Markovian, constitutes 
the major challenge in the general theory of quantum 
dissipation. ^^■^■^•^^ 

It has been shown that an exact reduced dynamics 
theory, in terms of hierarchically coupled EOM, does ex- 
ist in model Debye solvents that satisfies a semiclassical 
fluctuation-dissipation relationiSi2^ Based on this exact 
theory, we shall in this work construct an analytical rate 
expression for the simple ET system, without invoking 
such as resummation and timescale separation approxi- 
mations. As results, the present work will not just recover 
for kinetic rates the celebrated Kramers' turnoves^ii and 
Marcus' inversion^i^ behaviors, it will also reveal for ET 
thermodynamics such as Gibbs free-energy and entropy 
functions some interesting solvent dependent behaviors. 

The remainder of this paper is organized as follows. 
Section II treats an exact, nonperturbative theory of the 
reduced density matrix dynamics in Debye solvents. Af- 
ter a brief review of the hierarchical EOM formalism 
f Sec. lII At M^ we construct a continued fraction Green's 
fimction theory of quantum dissipation fSec. lIIBll . We 
further utilize it, together with Dyson equation tech- 
nique, to evaluate analytically the reduced dynamics of 
the simple ET system fSec. lII(J|) . Section Hill contributes 
to the development of reduced density matrix-based ET 
rate theory. Numerical studies in Sec. II VI will demon- 
strate not just for ET rates, but also for ET reaction 
(equilibrium) thermodynamics, their dependence on sol- 
vent environment. Finally, Sec.lVl concludes the paper. 



II. EXACT REDUCED DYNAMICS IN DEBYE 
SOLVENTS 

A. Hierarchical equations of motion formalism 

To describe the hierarchical EOM for reduced density 
matrix)2^i2^ let us recast the total ET Hamiltonian (eq^J 
in the stochastic bath interaction picture. 



HT{t) = H + H'{t), 



(3) 



with H and H'(t) representing the reduced system 
Hamiltonian and the stochastic system-bath coupling, re- 
spectively. 



H={E° + X)\b){b\ + V{\a){b\ + \b){a\), 



(4) 



(5) 



H'{t)^[U{t)-X]\b){b\. 

The stochastic solvation coordinate, 

U{t) = e'f'-*/^Ue~'^-^'^''e'^^'^''{hb - /ije-''*"*/'', (6) 

is assumed to be of Gaussian statistics. Thus, the effects 
of solvent on the ET system are completely determined 



by the solvent reorganization energy, 

A = (C/(0)^tr[t/(iK^] = ([/), (7) 

and the solvation correlation function, 

Cit-T) = {[Uit)-\][UiT)-\]). (8) 

In this work, we focus on the ET system in a Debye 
solvent (also called the Drude model), characterized by 
the following form of solvation response function^SS 

i{[U{t), C/(0)]) = -2ImC(t) ^ 2eit)hXje-'''. (9) 

Here, 8(t) is the Heaviside step function and 'y~^ = ti^ = 
T-D{eoo/sQ), with td being the Debye time parameter, 
and eo (^oo) the static (high-frequency) dielectric con- 
stant of the solvent. The corresponding solvation cor- 
relation function in the semiclassical high-temperature 
limit readsSS 



C(t) w A(2A:bT - ift7)e-T*. 



(10) 



For this model, the exact reduced density matrix dynam- 
ics has been constructed, in terms oS^iS^ 

Pn^ -[iL + n'-i)pn~iBpn+i-inApn-i\ ri>0, (11) 

which hierarchically couple the p = po oi primary interest 
and a set of auxiliary system operators {p„; n = 1, 2, • • • }. 
The initial conditions are Pn(0) = p(0)(5„o, and 



2XkBT 



(12a) 



AO 
BO 



.^7 



^2 Jfc)(6|,0]-z-^{|6)(6|,0}, (12b) 



[\b){b\,0]. 



(12c) 



Here, {•, •} denotes anticommutator. It has been shownSi 
that the individual auxiliary operator pn>o accounts for 
the 2n"^-order system-bath interaction contribution to 
the reduced dynamics of the primary interested p; see 
also the comments in the last paragraph of Sec. lHBI 



B. Continued fraction Green's function formalism 

Introduce the propagators {W„(t); n = 0, 1, • • • }: 
p„(t) = e-"^*Z^„(t)p(0); with UniO) = 5no- (13) 
Equations II II read now 

Unit) = -iUAn{t) - iBe-^'Un+i{t) - m^e^*Zi„_i (t) , 
which in the Laplace-domain are 

(s+iL)i^„(s) + i6Z^„+i(s+7)+my«^„_i(s-7) = (5„o-(14) 
Define the Green's functions {^^"•'(s);n > 0} via 

Uq{s) = g^'>Hs) = g{s), (15a) 

Ur,{s) = -ing^'^\s)AUn-i{s - -i)\ n>0. (15b) 



These equations will lead to 

1 



g^"Hs) 



s + iL + n(") {& 



n > 0, (16a) 



with 



n(") (s) = {n + l)Se("+i) (s + j)A. (16b) 

The above equations, which can be recast as 

71+ 1 



n(")(s) = S' 



s + 7 + 2L + n("+i)(s + 7) 



A (17) 



constitute the infinite continued fraction formalism for 
evaluating each individual n'^")(s) or t/("')(s). 

The Green's function t/('')(s) = t/(s) and its associated 
n(°)(s) = n(s) are of the primary interest. The former 
resolves the reduced density matrix evolution (cf. e allSai 
and c all 31 with n = 0), 

/"OO 

p{s)^ / dte-^'p{t)^g{s)p{0). (18) 

Jo 

This equation can be react as (cf. e all6al at n = 0) 



sp{.s) - p{0) = -iLpis) -n{s)p{.s), 



which in the time-domain reads 



(19) 



pit) = -iLp{t) - / dr n(i - t)p{t). (20) 

Jo 



Therefore, 



n(s) ==n(")(s) = 



dte-'* 



m, 



(21) 



represents the memory kernel in the Laplace domain. 

The initial input for the inverse recursive evaluation 
of n(s) fe allTII can be chosen based on the following 
observation. Each A is of second order in the system- 
bath coupling; thus the leading contribution of II'"' to 
the required 11 is of the {2nY^ order. Moreover, as the 
mathematical nature of continued fraction, convergency 
is also guaranteed practically for arbitrary strength and 
timescale of system-bath coupling. We can therefore set 
jj(Af4-i) _ Q^ with a sufficiently large N, to initiate the 
inverse recursive procedure, and evaluate n'-"'(s -I- nj); 
first at n = A^, then iV — 1, and so on, until the required 
n(o)(s) =n(s) is reached. 



C. Evaluation of tensor elements 

The tensor element of an superoperator (or Liouville- 
space operator) O is defined in the double-bracket nota- 
tion as^^i^ 



O 



jj' Mk' 



{{3j'\0\kk')), 



(22) 



so that 

(^= E 0,,.,kk'\jj')){{kk'\. (23) 

jj' ,kk' 

For a two-level system considered in this work, each ten- 
sor has 2^ = 16 elements. That O is Hermite implies 
Ojj'^kk' = 0*,j,,,^. Apparently, all n(") and ^(") are 
Hermite. 

To analyze the tensor elements of H^") fe all6bp . let us 
first examine A and B, defined by eas ll2bl and ll2cl respec- 
tively. They are found to be diagonal, with the nonzero 



elements of Aba ba 



-a: 



b.ab 



X{2kBT - ih-/)/h^, 



Abb.bb = -i2\^/h, and Bba,ba = -^Ib^ab = 1- As results, 
the only nonzero elements in II'^"^ ('e all6bp are 



(n) ^ TjCn) (n) ^ ttC") („) ^ -prCn) 






ba,bb' 



(24) 



and their Hermitian conjugate elements, and they are 
related to the Green's function tensor elements. 



^ — ^hn hnl -^ — ^hn nhT ^ — ^hn 1 



ba,ab^ 



'ba,bb^ 



(25) 



by [denoting 77 = A(2A:bT - ih'y)/h'^] 

x(")(s) = 77(7i-f l)x("+i)(s + 7), 

y(")(s) = -77*(n + l)y("+i)(.s + 7), 

z(")(s) = (7?-7?*)(n + l)Z("+i)(s + 7). (26c) 

To evaluate the involving Green's function elements 
via e all6al we apply the Dyson equation. 



(26a) 
(26b) 



g^g^g{iL' + ii')g, 



(27) 



with g being the diagonal contribution, and L and H' 
the off-diagonal parts of the involving L and 11 in e all6al 
respectively. Here and in eas l28cl and [551 follows, the 
common superscript (n) and argument s in both sides 
of equations arc implied. After some elementary algebra, 
we obtain 



X 



Y 



a^ 



■f3* 



|a + /?|2-|/?-y|2' 

|« + /3|2-|/3-y|2' 
1 



(28a) 
(28b) 



Z = — [{z-iV/h)X + {z* +iV/H)Y], (28c) 



with 



a = s + {i/h){E° + X)+x, (29a) 

(3 = s~^{Vlhf{2 + ihz/V). (29b) 

We have thus established from e all6al the expressions of 
{X, y, Z}(")(s) in terms of {a;, y, z}(")(s), which together 
with eas l26l constitute an analytical and exact formalism 
for the inverse recursive evaluation of the reduced dy- 
namics in Debye solvents. In the following section, we 
shall show that the ET reaction rate can be expressed in 
terms of {x, y, z}, i.e., the nonzero elements of dissipative 
memory kernel H in the Laplace domain. 



III. ELECTRON TRANSFER RATE: 
REDUCED-DENSITY-MATRIX FORMALISM 

We are now in the position to construct the reduced 
density matrix approach to ET rates. Let us start with 
p{s) [e alT^ . where p{t = 0) = \aa)). By separating p{s) 
into population vector P — [paa, Pbb]'^ ^^nd coherent vec- 
tor [pabjPba]'^ componcnts, and then using e all9l to ehm- 
inate the latter, we obtain the ET kinetic equations in 
Laplace domain as 



{i/h){E° + A) + x{s) is the same as 



with 



sP(s)-P(0) = A'(s)P(s), 



K{s) = Tpc(s + Tcc)"'Tcp - Tp 



(30) 



(31) 



Here, Tpc, Tec 2^cpj and Tpp denote the coherence- 
to-population, coherence-to-coherence, population-to- 
coherence, and population-to-population transfer matri- 
ces involved in e all9l respectively. Tensor analysis (cf. 
Sec. lII C|l results in Tpp ~ 0, and 



Tp, 



iV 



-1 1 
1 -1 



Ta 



T^or 



n 
n. 



10,11 



Tec =i{E° + A) 



-1 
1 



TT* TT* 

-"-•-lO.lO -'-'-10,01 

nio,oi nio,io 

Note that e al3()l in time domain reads 



P{t) = / dTk{t - t)P{t). 



(32a) 



(32b) 



(33) 



Thus K{s) is the resolution or the Laplace-transform of 
the memory rate kernel K{t). The total population con- 
servation implies the relation Kaj + Kj^j = 0; thus, c al33l 
is equivalent to 

Pa{t) =- f dTk{t-T)Pa{T)+ f dr k' {t ^ t) Pb {t) . (34) 

Jo Jo 

The forward and backward rate resolutions are 

k{s) = -Kaa{s)= / dte-'*k{t), (35a) 

"'0 
/•oc 

k'{s) = Kab{s)^ / dte-''k'{t). (35b) 



Together with e al31l easl32land Tpp = 0, we obtain 

fc(.) = S^Re. "(^) + ^(^) 



and 



k'{s) 



h? ^^|a(s)P-|y(s)|2' 



2\VW\a{s) + y{s)][l-ihz*{s)/V] 



h? 



-Re 



\a{s)\^-\y{sW 



(36a) 



(36b) 



a;(0), 



^ba,ab = y^^\ and 



Here a[s) = s 

e al29al while x = l-i-ba.ba = x"-"', y 

z = Hba.bb = -z^^-* are the same as easl24lat n = 0. 

The above formulations relate the ET reaction rate res- 
olutions k(s) and fc'(s), to the memory dissipation kernel 
resolution n(s); cf. easll9landl2()l The analytical expres- 
sions for the inverse recursive evaluation of the H(s) ten- 
sor elements had been detailed in Sec. lH Bl and Sec. lH CI 
We have thus completed the analytical expressions for 
the frequency-dependent rates of ET in Debye solvents. 



IV. RATES VERSUS THERMODYNAMICS: 
NUMERICAL DEMONSTRATIONS 

For numerical study, we focus on the ET rate con- 
stants, k = k{s — 0) and k' = k'{s = 0), which amount 
to the integrated memory rate kernel (cf. e al35|l . It is easy 
to show that in the weak transfer coupling (small V) and 
slow solvation (small 7) limit, the present rate constant k 
fe al36al at s = 0) recovers the celebrated Marcus' ET rate 
expression of eq|21 We shall also be interested in the re- 
action Gibbs free-energy, AG" = —kBT\n{k/k'), entropy 
A^" = -dAG°/dT, and enthalpy AH° = AG°+TAS°. 
We shall demonstrate how the rate constants and reac- 
tion thermodynamics functions depend on the model pa- 
rameters, reaction endothermicity E°, solvent reorgani- 
zation energy A, and longitudinal relation time tl = I/7. 
The other two parameters of the model are set to be 
T = 298 K and V = 1 kJ/mol, unless being further spec- 
ified. In fact, the value of V can be considered as the 
unit that scales other parameters used. The tempera- 
ture should also varies around T ~ 298 K in order to 
numerically evaluation of entropy. 

Note that in optical spectroscopy, one often uses a di- 
mensionless parameter, k~^ = Ti^^/Tk^TX/h, to measure 
the slow [k, ^ 1) and the fast {n ^ 1) solvent modula- 
tion hmit. For A = 3 kJ/mol at T = 298 K, that k, = I 
corresponds to tl — 16.5 fs, while the typical ET sol- 
vation correlation timescale is of picoseconds. We will 
show that the parameter k can also be used in ET rate 
problem. 

Figure Q] depicts ET rate k as the function of solvent 
relaxation time tl ~ I/7 at four representing values of 
endothermicity, E° = 0,-1,-3, and —5 kJ/mol. The 
solvent reorganization energy A = 3 kJ/mol. In general, 
the ET rate consists of the barrier crossing and the coher- 
ent tunneling contributions. When E° = (symmetric 
case), the system is in Fermi resonance and the ET is 
dominant by coherent tunneling. The observed rate in 
this case exhibits the motional narrowing behaviorsiiS^ 
the faster the solvent modulation is, the larger the co- 
herent resonant tunneling rate will be. In a nonsymmet- 
ric {E° ^ 0) system, the barrier crossing is significant 
and the observed ET rate exhibits clearly the Kramers' 
turnover behaviors^*^ 

The rates observed in the slow solvent modulation 
(large tl which amounts to large viscosity) region of 



Fig.n are closely related to the Marcus' inversion be- 
haviors. Figure |21 depicts the ET rate, in term of 
ln(fc/fci„ax), as the function of reaction endothermicity 
E°, in the slow modulation region (tl — 10 ps), at two 
specified values of transfer coupling strength, V = 0.01 
and 1 kJ/mol. When the transfer coupling is small 
{V — 0.01 kJ/mol), the rate in the slow solvent mod- 
ulation regime does have the Marcus' nonadiabatic form, 
ln(fc/fcmax) = -'iE° + A)V(4/cBrA) of eq[2 cf. the thin- 
solid vs. thin-dash curves. In the case of V^ = 1 kJ/mol, 
however, the rate deviates significantly from the Marcus' 
expression. The poorly-fitted parabolic function there, 
-iE°-E^^^)^/{C4:kBTX) (dash-curve), is found to be of 
-E^max = -2-4 kJ/mol and C = 0.3. 

Figure |21 shows the 3d-plot of reaction Gibbs free- 
energy AG° as the function of (A, tl), exemplified at 
T = 298 K with E° = -3 kJ/mol and F = 1 kJ/mol. Re- 
ported in Fig.01 are some representing 2d-slices of Fig.|3 
for AG°, together with the numerically evaluated reac- 
tion entropy AS° and enthalpy AH° . The observed fea- 
tures here are listed as follows. 

(i) The basic symmetry requirements, such as 
AG°{-E°) = -AG°{E°) which implies also AG°{E° = 
0) = 0, hold in general; 

(ii) AG°, AH° and E° are of same sign and \AH°\ > 
I AG" I > \E°\, implying that the enthalpy and the en- 
tropy play opposite roles on the reaction Gibbs free- 
energy; 

(iii) The dependence of AG° on (A, tl), as shown in 
Fig-El is qualitatively similar to that in the weak transfer 
coupling regime; 

(iv) In general, AG° approaches to a constant in both 
the fast (k ^ 1) and slow (k <C 1) modulation regimes. 
The smaller V is, the closer AG° to E° will be. For the 
system demonstrated in Fig.|21(or Fig.QJ where E° ~ —3 
kJ/mol with V = 1 kJ/mol, |AG°(k > 1)| > |AG°(k < 
1)1; however the sign could be opposite for small \E°\ 
systems (not shown here); 

(v) In the modest modulation regime {k ^ 1), the re- 
action thermodynamic functions (AG°, AH°, and A5°) 
exhibit certain nonlinear dependence on (A,rL). In par- 
ticular, the magnitudes of these ET reaction thermody- 
namic functions are of maximum values around k ^ 1 
when A> \E°\. 

Interestingly, the observed Kramers' turnover of rates 
occurs also around the k '^ 1 (intermediate friction) re- 
gion. The relevant forward and backward rate constants 



are given in Fig.|S| This may at least partially account for 
the nonmonotonic dependence of thermodynamics func- 
tions on the solvent environment in this region, as they 
are related to the rates via the detailed balance relation. 



V. SUMMARY 

In summary, we have constructed a formally exact, 
nonperturbative ET rate theory in terms of the resolution 
of the memory dissipation kernel (cf. eas l36l with eas l21l 
and l24|l . For Debye solvents in which the solvation corre- 
lation is characterized by an exponential function fe allO|l , 
the ET rates, or rate resolutions in general fcas l36|) . can 
be evaluated readily via the inverse-recursion formalism 
(cf. easl26lwith eas l28|l . Not only does it recover the cel- 
ebrated Marcus' inversion and Kramers' turnover behav- 
iors of the ET rates, the present formalism also provides a 
microscopic theory for the effects of solvent environment 
on the ET thermodynamics. The dependence of reac- 
tion thermodynamics on solvent environment is found 
to be quite dramatic, especially in the region where the 
Kramers' turnover occur(s). This observation suggests 
the possibility of utilizing the thermodynamics data to 
extract such as the solvation correlation time parameter 
(cf. the left panel of Fig.^l). 

This work has also developed a nonperturbative the- 
ory of reduced density matrix dynamics, in terms of con- 
tinued fraction Green's function fSec. lIIBll . This for- 
malism can be used together with Dyson equation tech- 
nique for efficient and analytical evaluation of reduced 
dynamics fSec. lIICJl . The present formalism is exact for 
the Debye solvents fe allO() . However, as the semiclassi- 
cal fluctuation-dissipation theorem is involved in c all 01 
the reduced density matrix and rates may become nega- 
tive if fceT < (|1/|2 + i|£;°|2)i/2. Generalization of the 
present continued fraction Green's function approach to 
non-Debye solvents at arbitrary temperature is feasible 
and will be developed in future. 
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FIG. 1: Electron transfer rates as the functions of solvent 
longitudinal relaxation time tl = I/7, for some specified 
values of endothermicity E° . The solvent reorganization en- 
ergy, transfer coupling strength, and temperature are A = 3 
kJ/mol, V = 1 kJ/mol, and T = 298 K, respectively. 



FIG. 2: The normalized rates, in the slow modulation regime 
(tl = 10 ps), as the functions of E° . The Marcus' parabolic 
relation, ln(fc/fcmax) = -{E° + Xf/i-iksTX), where A = 3 
kJ/mol, is recovered in the case of V = 0.01 kJ/mol, but not 
for V — 1 kj/mol. The latter case neither fits well with a 
parabolic function (dash curve); see the text for details. 



FIG. 3: Reaction Gibbs free-energy AG° as the function of 
solvent parameters (A,tl), for an ET system with E° = —3 
kJ/mol and 1/ = 1 kJ/mol at T = 298 K. 



FIG. 4: Reaction thermodynamics functions AG", AS° , and 
AH°: (a) as the functions of A at some selected values of tl; 
(b) as the functions of tl at some selected values of A. The 
ET system is same as that of Fig.|21 



FIG. 5: The forward and backward rate constants, k (upper 
panels) and k' (lower panels), relevant to the reaction ther- 
modynamics functions in Fig.^] 
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